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To investigate correlation effects on antiferromagnetic order in Fe pnictides, we apply a variational 
Monte Carlo method to a two-orbital model. We obtain a small ordered moment consistent with experi- 
mental observations even for a Coulomb interaction comparable to the band width. Studies of estimation 
of the Coulomb interaction for Fe pnictides suggest values comparable to or slightly smaller than the band 
width, and much larger ordered moments have been obtained by the Hartree-Fock approximation for such a 
large Coulomb interaction. Thus, the correlation effect is important for Fe pnictides at least quantitatively. 
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The discovery of superconductivity in LaFeAsOi-^F^ 
with a high transition temperature T c = 26 K 1 ' has stimulated 
extensive and intensive studies on Fe pnictides. Superconduc- 
tivity takes place around the magnetic phase boundaries '~ 5 ' 
as in high-T c cuprates. Such a similarity suggests that mag- 
netism is playing an important role in the emergence of su- 
perconductivity, and it is highly desirable to unveil the micro- 
scopic origin of magnetism characteristic to Fe pnictides. 

To unveil the magnetism in Fe pnictides, the present au- 
thors applied Hartree-Fock approximation 6 ' to a two-orbital 
model. 7 ' The results are summarized as follows. The antifer- 
romagnetic order with ordering vector (tt, 0), in the unfolded 
Brillouin zone with one Fe ion per unit cell, is stabilized by 
the nesting between hole and electron pockets. This antifer- 
romagnetic state inevitably accompanies ferro-orbital order, 
since the ordering with (tt, 0) breaks the equivalence of x and 
y directions, and as a result, the occupancies of d zx and d yz 
orbitals become different. Under such ferro-orbital order, the 
lattice should be distorted from a tetragonal to orthorhombic 
structure through an electron-lattice interaction. Even in the 
antiferromagnetic state, a band gap does not open at some 
points in the Brillouin zone due to multiorbital nature of the 
bands. Therefore the system remains metallic in the ordered 
state. These results are consistent with experimental obser- 
vations of magnetic order with (tt, 0), lattice distortion, and 
metallic conductivity. 

However, in our previous Hartree-Fock result, the ordered 
moment is large in contradiction with experimental observa- 
tions, 5, 8 ~ 14 ' if we take a large Coulomb interaction compara- 
ble to the band width. Studies of estimation of the Coulomb 
interaction for Fe pnictides suggest values comparable to 15 ' 
or slightly smaller than 16 ' the band width. For such large val- 
ues of Coulomb interactions correlation effects beyond the 
Hartree-Fock approximation may be important. In particular, 
correlation effects are expected to reduce the magnitude of the 
ordered moment. Indeed, importance of the correlation effects 
is discussed for a three-orbital model by using a Gutzwiller 
approximation. 17 ' 

In this paper, we investigate correlation effects on mag- 
netism by applying a variational Monte Carlo (VMC) method 
to the two-orbital model. While the VMC method has been 



applied to a five-orbital model with a partially-projected 
Gutzwiller wavefunction, 18 ' only Fermi-surface distortion 
and superconductivity are discussed there. In the VMC 
method, we consider a Gutz wilier-projected wavefunction as 
a variational wavefunction. We show that this wavefunction 
contains substantial correlation effects beyond the Hartree- 
Fock approximation while this wavefunction is simple enough 
for numerical calculation. 

In the two-orbital model, we consider a square lattice of Fe 
ions with d zx and d yz orbitals. 7 ' 19 ' The model Hamiltonian is 
given by 

h.r,r',a i,r 

where Ci Ta is the annihilation operator of the electron at site 
i with orbital r and spin a (=t or I) and Cfc T(T is the Fourier 
transform of c; TCT . The orbital indices t = x and y repre- 
sent d zx and d yz orbitals, respectively. The number operators 
are defined by n lTa = c\ T<J c iTIJ and n lT = X) CT n iT0 .. The 
coupling constants U, U', J, and J' denote the intraorbital 
Coulomb, interorbital Coulomb, exchange, and pair-hopping 
interactions, respectively. The relations U = U' + J + J' 
and J = J' hold for the t2 g orbitals, 20 ' and we use them. 
We use the hopping parameters proposed by Raghu et a/. 7 ' 
and the coefficients in the kinetic energy terms are given 
by tkxx — —2ticosk x — 2i2Cos/cj, — 4t3Cosk x cosk y , 
£ fcyy = —2t2COsk x — 1t\ cos k y — 4^3 cos k x cos ky, and 
£kxy = £kyx = — 4^4 sin k x sin k y , where t\ — —t,t2 = 1.3t, 
t% = t4 = — 0.85t, and we have set the lattice constant unity. 
The band width is W = 12i. 

We consider the variational wave function given by 

I*) = Pcm = nt 1 - (! - ff7)l*r>H]|*>, (2) 

where Pc is the Gutzwiller projection operator for onsite 
density correlation. 21_24 ' denotes projection onto the 
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state 7 at site i and g 1 is the variational parameter controlling 
the probability of state 7. There are sixteen states at each site 
in the present two-orbital model. The Hartree-Fock type wave 
function |$), which describes a charge, spin, orbital, and spin- 
orbital coupled ordered state, is given by 



i*> = n 6 £ } >>' 



(3) 



where a is a band index and |0) is the vacuum. The quasi- 
particles occupy N a states for each spin a from the lowest 
quasiparticle energy state, where N a is the number of elec- 
trons with spin a. Here we consider a half-filled case, and we 
set JVf = 2V| = N, where N is the number of the lattice 
sites. The quasiparticle states are obtained by diagonalizing 
the following 4x4 matrix: 
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where Q = (tt, 0) is the ordering vector. The quasiparticle 
gap in the ordered state is given by 



A Ta = A (S rx - S Ty ) + A so (6v t - <5 CT 4.)((5 r2 



Sry), 
S T y) 



(5) 



(6) 



A t<j q =A c q + A sQ (S a1 - - 8 a i) + A q(* 

+ A soQ (S a ^ - 6 a _i)(5 TX - 5 Ty ), 

where A and A so denote the gaps for uniform orbital and 
spin-orbital ordered states, respectively. A c q, A s q, A q, and 
A so q denote the gaps for antiferro-ordered states of charge, 
spin, orbital, and spin-orbital, respectively. We also take them 
as variational parameters. 

For this variational wavefunction, we evaluate energy by 
the Monte Carlo method, and optimize variational parameters 
to find the state which has the lowest energy. We set all A ra 
and A Ta Q zero to evaluate energy of the paramagnetic state, 
that is, we optimize only Gutzwiller parameters g-y. For the 
antiferromagnetic state, we also vary A D , A s q, and A so q. We 
also evaluated energy by varying all A T(J and A T(7 q for some 
values of U, but we could not find a solution which has lower 
energy than the antiferromagnetic state. The calculations are 
done for an 8 x 8 lattice with an antiperiodic boundary condi- 
tion for both directions. 

Figure 1 shows energy as functions of U obtained with the 
Hartree-Fock approximation 6 ' and the present VMC method. 
The energy is lowered by the correlation effects beyond the 
Hartree-Fock approximation. 

Figure 2 shows the ground state energy as a function of U 
measured from that of the paramagnetic state. The transition 
from the paramagnetic state to the antiferromagnetic state oc- 
curs at U > It. If the energy difference is proportional to 
U 2 around the transition it is of second order, and if the en- 
ergy difference is proportional to U the transition is first order. 
However, it is difficult to distinguish a second order transi- 
tion from a weak first order transition as is obtained by the 



J=0.1U 



Hartree-Fock 
VMC: 8x8 



6 

13 It 



8 



10 12 



Fig. 1. (Color online) Energy E as functions of the Coulomb interaction 
with J = 0.1(7 obtained with the Hartree-Fock approximation 6 ' and the 
VMC method. 
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Fig. 2. (Color online) Energy Eaf of the antiferromagnetic ground state 
measured from energy Ep ar n of the paramagnetic state as a function of U 
with J = 0.1(7. 



Hartree-Fock approximation 6 ' from the present results due to 
numerical accuracy. 

Figure 3 shows the ordered magnetic moment m S Q evalu- 
ated for the optimized wavefunction. m S Q is defined as 



(7) 



where denotes the position of site i and (• • • } represents 
the expectation value. To check the finite size effect of the 
model, we also show the results for a 10 x 10 lattice. The 
finite size effect on m S Q is weak in particular for the large 
m S Q region. In the result of the Hartree-Fock approximation, 
there is a small but finite jump in to s q, 6 ' while it is invisible in 
the scale of Fig. 3. For the results of the VMC, as in the energy 
difference, it is difficult to determine whether the transition is 
first order or second order. If it is a first order transition, the 
jump in the magnetic moment at the transition is very small. 
The ordered moment is not large for U < 9t. By comparing 
the results by the Hartree-Fock approximation and the present 
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Fig. 3. (Color online) Ordered moment m s Q in the antiferromagnetic 
ground state as functions of U with J = 0.1U obtained with the Hartree- 
Fock approximation 6 ' and the VMC method. Note that if a fully polarized 
state is realized, it has m s Q = 2. 

VMC results, we conclude that the correlation effect strongly 
reduces the value of the ordered moment and such an effect is 
important for Fe pnictides. 

Around U = 9.6i, we find another phase transition within 
the antiferromagnetic phase. This phase transition is of first 
order and it is probably a metal to insulator transition, since 
the energy gain by the kinetic energy is reduced at U > 9.6i 
(not shown). 

We have also searched for a ferro-orbital ordered state with- 
out antiferromagnetic order, that is, the gap parameters are set 
zero except for A Q , but we could not find such a state as a 
ground state. In the antiferromagnetic state, the order param- 
eter Too = (1 /AT) J2i a ( n ix(T — niya) for the ferro-orbital order 
becomes also finite due to symmetry lowering. However, the 
values are too small and we cannot determine to q confidently 
due to our numerical accuracy. 

To summarize, we have applied the variational Monte Carlo 
method to a two-orbital model to investigate correlation ef- 
fects. Then, we have found that the ordered moment in the 
antiferromagnetic state is strongly suppressed by the corre- 
lation effect. Thus, to obtain a small ordered moment as in 
experimental observations, for U < W, we should take cor- 



relation effect into account properly. 
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